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ABSTRACT 


New genetic operators are described that assure preservation of the feasibility of candidate solutions to 
any discrete and linearly constrained optimisation problem. The design of these operators is the result of 
extensive theoretical investigations, with particular assiduity devoted to considering the most challenging 
examples of this type. 

Attention is largely centred on problems that have defied satisfactory solution by traditional means, 
because of poorly behaved or imprecise objective functions, high dimensionality, or intractable algorithmic 
complexity. Evolutionary Algorithms, inspired by adaptation in biological systems, are ideally suited to 
such arduous conditions, as testified by their increasing popularity and successful exploitation in 
resolving numerous difficult problems. 

Their obvious attraction notwithstanding, the applicability of evolutionary algorithms has suffered by the 
deficiency of general techniques to manage constraints, a feature common to many problems, and the 
asperity is compounded when some proportion of the decision variables are discrete. The enhanced 
operators presented here guarantee the feasibility of proffered aspirant solutions with respect to a system 
of linear constraints. 

Exploration of highly constrained systems indicates the possibility of performance degradation, with the 
diminishing probability of discovering a feasible and meaningful exchange of information between 
candidate solutions. Relaxation of the crossover operator is proposed to alleviate this, and the effective 
utilisation of this in the context of the overall algorithm suggests a modification in which the population is 
permitted to transiently increase in size. 
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Feasibility and Genetic Algorithms: 

The Behaviour of Crossover and Mutation 

EXECUTIVE SUMMARY 

This document details the results of enabling research aimed at providing advanced 
techniques for building sophisticated automated planning and decision support tools. 
In particular, a broad range of apparently dissimilar resource management problems 
such as surveillance planning, resupply distribution planning, and fire support 
planning are, in fact, intimately related. 

Solving these kinds of problems can be deceptively difficult; they may require 
unrealisable computational resources to resolve problems of even modest size. In 
addition, some of these planning and decision support functions display complicated 
interrelationships between the factors that together define the relative desirability of 
one possible solution over another. This nonlinearity itself poses substantial barriers to 
constructing effective tools for improving the timeliness and quality of decisions. 

Additional considerations arise from the nature of the context of military operations 
and the translation of the military problem into mathematical terms; some subtle and 
poorly understood factors may not be fully represented in the mathematical 
formulation. The accepted resolution to this is to augment the ability of the human 
decision-maker with automated assistance, rather than to attempt to completely 
replace the human with a fully automated tool, at least in the first instance. To achieve 
this, the proposed solution technologies should provide the human with a range of 
plausible solutions to consider before he or she implements a final decision. 

In response to these realisations, research has focussed on discrete optimisation 
problems and the discrete component of mixed optimisation problems that defy 
satisfactory elucidation by more traditional means. Attention is particularly focussed 
on problems of high dimensionality that describe poorly behaved objective functions, 
and those for which fully deterministic algorithms present intractable algorithmic 
complexity. It is also presumed that an appropriate method of solution should develop 
a range of qualitatively comparable candidate solutions close to the theoretical optima. 

Evolutionary and genetic algorithms are a class of randomised search methods that are 
squarely aimed at problems of this character. Inspired by the adaptive behaviour of 
biological organisms, all evolutionary and genetic algorithms maintain a population of 
potential solutions, with each new generation developed from the last using a set of 
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stochastic transition operators. This contrasts sharply with traditional approaches, 
which typically produce a sequence of single solutions deterministically. 

Despite their inherent robustness and obvious utility, the application of evolutionary 
and genetic algorithms to resource management problems suffers from their inability 
to adequately manage systems of constraints. This deficiency is compounded further 
with the need to constrain decision variables to integral values. This document firstly 
examines the behaviour of usual genetic crossover and mutation operators with respect 
to any system of constraints, by describing the probabilities that the generated 
solutions lie within the feasible region that such a constraint system defines. Building 
further upon this basis when the constraints are all linear, enhanced operators are 
proposed that guarantee the feasibility of their results. 

Highly constrained problems are characterised as those that possess a very small 
feasible region with respect to its smallest encompassing interval. Investigation of the 
behaviour of genetic operators in such a setting reveals the potential for greatly 
degraded performance as the probability of generating only feasible candidate 
solutions collapses. However, the probability that at least some of the solution 
proffered by crossover and mutation may remain significant. This suggests a scheme in 
which the population is permitted to fluctuate in size. 
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1. Introduction 

The Evolutionary and Genetic Algorithms are a class of stochastic relaxation 
techniques, based on the progressive adaptation of biological systems, that are 
attracting increasing interest for the resolution of a wide variety of challenging 
applications [1-12]. In contrast to sequential search methods, a genetic algorithm 
maintains a large population of candidate solutions, rather than a single iterate, and 
each population is developed by the application of stochastic transition operators to its 
predecessor. 

Genetic algorithms also differ in the character of the development of new points for 
exploration; while most traditional methods act directly on the value of an iterate, the 
genetic operators utilise an encoding of the solution, in the form of a string in some 
appropriate alphabet. This string is the 'chromosome' that defines an individual 
population member, as distinct from its actual value, or 'phenotype'. 

These operators can be grouped according to their biologically-inspired functions: 
selection operators define the virility with which each population member reproduces, 
crossover operators describe the mechanism by which individuals mate to produce 
offspring, and mutation allows an occasional alteration in genetic structure. 
Historically, the term 'Evolutionary Algorithm' has related to methods that utilise 
selection and mutation alone, while 'Genetic Algorithms' also employ crossover, and 
sometimes a variety of other additional mechanisms. Here, the distinction will not be 
vigorously supported, as the discussion presented permits either. 

The attraction of the genetic and evolutionary algorithms is enhanced by their 
conceptual simplicity, and the capacious account of disparate problems that have 
yielded to this kind of approach. Although its exploitation to the enhancement of 
overall performance is possible, no gradient or other additional information about the 
objective function need be available, so that noisy, stiff, and non-differentiable 
problems in principle pose no additional challenge. Furthermore, the maintenance of a 
large population of candidate solutions reduces the tendency to converge to local, but 
not global, optima, rendering this approach a vigorous contender in the illation of 
multimodal problems. 

Despite their obvious allure, the applicability of genetic algorithms has suffered in the 
lack of adequate and sufficiently general techniques for the management of systems of 
constraints [3,6-11]. In most existing approaches, the genetic operators are permitted to 
generate infeasible candidate solutions, and these are penalised by modifications of the 
objective function (penalty functions), or later corrected (repair algorithms). The 
adoption of highly specialised representations of the problem (decoder algorithms), 
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that ensure, or magnify the probability of, the manufacture of feasible solutions, have 
also been suggested. 

Penalty function methods [3,9-10] discriminate against a possible solution according to 
its violation of the constraint system. Regardless of the magnitude of these penalties, 
the result may be calamitous, particularly if the number of infeasible solutions 
generated in each iteration is large. Both decoder and repair algorithms are frequently 
discovered to be computationally intensive, difficult to design, and highly problem- 
specific [6,11]. Furthermore, a decoder that still permits the possibility of developing an 
infeasible solution necessitates the utilisation also of penalties or repair algorithms, 
enhancing further the undesirability of this approach. 

In response to these realisations, awareness has more recently been focused [6-8] on the 
alteration of the genetic operators themselves to preserve feasibility. In doing so, the 
potential of the constraint system to reduce the size of the search space is fully realised, 
facilitating the possibility of greatly improved behaviour. The additional cost imposed 
by these modifications is more than offset by the would-be burden of penalty 
functions, decoders, or repair algorithms. 

It is a new methodology, conforming to this philosophy, for the resolution of problems 
constrained to integer values and by a system of linear inequalities that is described 
here. However, it should be understood that this might be easily combined with 
methods for continuous variables presented elsewhere [6], for the solution of mixed- 
integer type models. Within this context, then, the aspiration is the minimisation or 
maximisation of some function 


/ : Z" -> R, 


subject to the constraint system 


A- x <b , l < x < u , jceZ", 

where A = [a i j ]e R mx " and b = [b t ] e R m , with l = [l } ] e Z" and u = [uj]e Z" . 

Note also that no equality constraints are included; these may be converted to a system 
of this form with a reduction in the number of variables, and therefore also in the size 
of the search space [6-8]. 


2. Two-Point Crossover 


2 


The genetic selection operator is the emblematisation of the process of natural selection 
in biological systems, in which individuals survive to reproduce at varying rates [1,3- 
5,7-8,11]. The genetic algorithm emulates this discrimination by constructing a 'mating 


DSTO-TN-0320 


pool’, with each individual represented according to its 'fitness', which is its objective 
function value compared in some sense with those of the rest of the population. Pairs 
of individuals will then be chosen at random from this mating pool to undergo 
crossover. 

Modelled on the reproduction of biological organisms by the exchange of genetic 
information, the crossover operator facilitates the parturition of new individuals from 
old. Crossover speculates on new solutions constructed by combining fragments of old 
solutions that have previously worked well, according to their fitness values. The 
emergent population is therefore an innovation comprised of solutions that are the 
juxtaposition of fragments that have been successful in the past. 

The two-point crossover achieves this by swapping the components of two vectors in 
some random interval. For each pair of parent solutions chosen from the mating pool, 
two randomly generated indices mark the boundaries of the corresponding segments 
of information to be inter-changed. 

In seeking to understand the ramifications of restricting this operator to the 
development of feasible solutions, consider the feasible region, which is a convex 
polyhedron P in n-dimensional space, 

P = {teZ" . A - x <b a l < x < u\. 

This set is covered by the smallest bounded set AcZ" that forms a closed interval 
around P , 


AEA^r^eZ" :I < x <u\, 


where 


/ = [Jj • /,■ = mini*, : x = [x j ]" =1 AAx<bAl<x<u\, 


and 


u = [«,•]?= i • «, = max{c ; : x = [x ]" =) a A- x < b a l < x < u\, 

X J J 

are the vectors whose components are the smallest and largest, respectively, 
corresponding components of all feasible solutions x. 


The crossover operator takes as its arguments two vectors in n-dimensional space and 
the indices defining the crossover segment, and produces from this a new pair of 
vectors. No order is defined between the vectors of such a pair, and therefore all such 
pairs is a set 

n = ix,y>:x,yeZ"}. 
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The indices defining the information exchange indicate some number of consecutive 
corresponding components of the argument vectors; the set of all possible indices is 
described by the set 


X=§:l<i<nAie N}. 

Two-point crossover is therefore a function 

E :ilxx x X & 


defined as 


f{x',y'} if h < k 
E(Uy},M)= \ 7 . 

[ {x, y) otherwise 

and where x, y , and the resultant vectors x' , y' have the form 

x = [*\ ■•■x n f, y = \.y\---y n ] T > 

x’=[x r --x h y h+ r--y k x k+] ■■■x n ] T , 
y-[yi —y* x h+i —x k y k+ r-y n f ■ 

The crossover is called 'meaningful' if h < k ; when h > k no exchange of information 
occurs. Furthermore, the set of application of crossover will, for convenience, be 
abbreviated as follows: 


E lu (x,y) = E({x,y},h,k). 

Any meaningful crossover for which E h k (x, y) = {x, y} will be termed 'trivial', including 
the case when x = y. All other meaningful crossovers, having the potential of 
generating points previously unseen, are 'non-trivial'. 

Clearly, the crossover of feasible solutions cannot result in the production of points 
outside the surrounding interval AcZ” (note here also that the converse is not true). 
Therefore, in considering the result of a crossover operation, attention may be entirely 
restricted to this interval. In the presumption that the points forming the old 
population are all feasible, the crossover will only be applied to elements of the feasible 
region P . That is, only crossovers of the form 

E h,k (x,y)» ({*, y},h,k) e Of/ 5 ) 
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are relevant, where O(P) is the set of all meaningful crossovers between points in any 
set 5 such that Sc A, 

0(5) = ^[x,y),h,k): x,y e S a x ^ y Ah,k e x < k}. 

The number of meaningful crossovers between the points of a set Sc A is easily 
derived [7]. 


Lemma 2.1 [7]: 

The number of meaningful crossovers between solutions in 5 c A is 

10(5) I = — (« 2 - «)• 15 l (l5 l-l). 

4 

In particular, when 5 = P , this yields the number of meaningful crossovers possible 
between feasible points. 

The potential exists for the extemporaneous generation of infeasible candidate 
solutions, whenever the number of points I A I in the enclosing interval A transcends 
the number I P I in the feasible region P . The design of an appropriate mechanism for 
managing the crossover depends upon exactly how great this propensity is. 

The most immediate strategy is to permit only those crossovers that would 
manufacture nothing but feasible points; crossovers of this type, namely 

\({x,y},h,k)e O(P) AE hk (x,y) c P\, 

will be termed 'feasible', or 'entirely feasible' in emphasis of the fact that only feasible 
solutions are manufactured. 

The number of crossovers that produce only feasible points, from points in P, 
potentially bears considerable impact on the overall performance of the genetic 
algorithm. If the likelihood of a successful search for a crossover segment that 
preserves feasibility is remote, then little progress towards convergence is possible. 


Theorem 2.2 [7]: 


The probability that a crossover randomly chosen from <t>(P) is entirely feasible is 




I PI-(IP 1-1) 

I A I -(I A I —1) " 
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A graph of this probability distribution appears in [8]. Put another way, of the 
■j{n 2 - n)- 1 P I (I P I -1) meaningful crossovers between feasible points, 

M(P)=Pr(s w (^y)c/>)lO(P)l 

1, 2 x I PI 2 (IPI-1) 2 (!) 

= — (rt -n) - 

4 I A I (I A I -1) 


of them are entirely feasible. 

As the size of the feasible region P diminishes with respect to that of the enclosing 
interval A, so too does the probability (Theorem 2.2) of a crossover being feasible. The 
result is a decrease in the number of entirely feasible crossover operations (see 
equation (1)), and consequently in the performance of the genetic algorithm overall. 

The amelioration of this is facilitated in the realisation that not all crossovers producing 
feasible solutions are themselves entirely feasible. The class of crossovers that produce 
only feasible solutions is contained within the wider set of crossovers that produce at 
least one feasible point. By exploiting this, at least one new solution might be obtained, 
where none could be illuminated before. 

The set of 'half-feasible' crossovers, that produce only one feasible point, is 
i{x,y},h,k) e O(P) :l S M (*, y) fl P 1= l}, 

and its size can be ascertained similarly to that of the entirely feasible crossovers, by 
first considering the probability that a crossover between points in P belongs to this 
set. 


Theorem 2.3 [7]: 

The probability that a crossover randomly chosen from ®(P) is half-feasible is 

Pr( I t (,.y) fl f I = l)= ———— ■ 
v ' ' IA1-0 A I—I) 

A graph of this probability distribution can be found in [8]. The number of half-feasible 
crossover operations, with respect to the feasible region P, is therefore 


H(P)=W^ CP'-OCA'-'P 1 ), 

2 IA I (I A I —1) 


( 2 ) 
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Comparing Equations (1) and (2), the number of half-feasible crossovers is discovered 
to surpass that of the entirely feasible variety, when the feasible region consumes less 
than about two-thirds of its enveloping interval: 


_ — 2 
H(P) > M(P) whenever 1 P k — (I AI +1). 


(3) 


Under such conditions, the utilisation of the half-feasible crossovers may provide 
considerable advantages over restricting sedulity to those that produce only feasible 
solutions. However, tire realisation of this potential to the greatest advantage itself 
demands careful application. 


3. Maintaining Feasibility 

In order to afford the limitation of the choice of crossover segment to those that 
produce only feasible solutions, a sensible and efficient search is required. The 
exploration of all potential crossovers, for a given pair of points, is prohibitively 
expensive; first appointing one of the segment boundaries, by a random selection, 
reduces the number of possibilities from 0(n 2 ) to O(n). 

The actual search itself might then be effectuated by producing, in turn, each of the 
would-be resultant individuals, and appraising their satisfaction of the constraint 
system. However, in the case of a linear system, as prescribed, this disregards the 
relationship that exists between crossovers that have segments with one boundary in 
common. More thorough excogitation reveals that the entire search for an entirely 
feasible crossover between two given vectors can be made in a total of 0{mn + n 2 ) time. 

For brevity of notation, the residual r : Z" —> R m of a vector in n -dimensional integer 
space is a function defined as 


r(x ) - b - A- x , 


representing the amount by which a solution vector x violates or satisfies the 
constraints. The proposal is made that the residual of each vector be stored with it, to 
enable a rapid search for feasible crossovers; the residuals of the new vectors produced 
can be quickly found from those of the parents. The search method itself is justified by 
the following theorem. 


7 
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Theorem 3.1 [7]: 

Consider the crossover {*',/} = E hk (x,y), with x' and y' as described in the definition 
of crossover presented in Section 2. Then 

x’e P <=> 8 (h ' k) < r(x) , and y'e P <^>-r(x) <8 (h ’ k) , 
where 8 (h,k) = 8 (h ' k) (x, y) satisfies 

5 ^=8 (h ’ k) + (y h -x h )-[a ith ]l u 
S (hMl) = <5 (W) + (y k+ 1 ~ x k+l ) • [a iMl ]“,, 


$<*■*> = 0 , 

for any h,k:l<h<k<n. Furthermore, 

r(x') = r(x)-8 (h,k) and r(y') = r(y) + S (h ' k) . 

Given a pair of vectors from the mating pool, the quest for a feasible crossover begins 
with the random choice of the first segment boundary. This may be interpreted either 
as the lower boundary h , or as the upper k , so the search proceeds in two stages. All 
feasible crossovers that have this chosen point as the lower delimitation h might be 
generated by initialising 


S (h ’ k) =0 

and successively applying the recurrence relation 

8 khMl) = * (W) + (y k+l - x k+] ) • [a iMi ]£,, 

with each such crossover segment producing only feasible candidate solutions if 

-r(y)<8 (h ’ k) < r(x). 

A small example illustrating this procedure with a system of three linear constraints in 
five dimensions appears in [8]. 

The search, applying instead the alternative recurrence stated in the theorem, for 
crossovers with the chosen point used as the upper boundary instead, may be then 
performed. Early termination of the whole procedure can be accomplished by 
probabilistically deciding to accept or reject each feasible crossover, as they are 
illuminated. Given that there are a total of n -1 possible choices of crossover segments 
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that have any particular index as one of the boundaries, the probability that any 
random crossover in O(P) is entirely feasible (Theorem 2.2), and presuming a uniform 
distribution in the choice of the crossover from all those discovered to be feasible, the 
probability (on average) that a feasible crossover is chosen and applied [7] is 

__ 1 _ 

Pp (n-l)-Pr(E w (r,}) c P) 

This assumes also that the expected number of crossovers found is at least one; when 
this condition is violated, the feasible region is exceedingly small, demanding instead 
the exploitation of a modified scheme, to be discussed later. 

In most practical situations, the parameter p F cannot be known prior to the 
application of the genetic algorithm, as it is defined in terms of the size of the feasible 
set. By keeping a count of the total number of feasible crossovers found t F , and the 
total number of possibilities explored t c , however, p F can be approximated using 

Prfc/a (*>>’) c:P)~y-- 

l c 

The estimation of p F used in the algorithm is periodically updated, typically at the 
end of each generation (major iteration). As the algorithm progresses, the greater the 
number of crossovers it samples, so that this approximation rapidly converges to the 
actual value. 

Given the dominance of crossovers that manufacture a single feasible point over those 
that produce two, in more highly constrained examples (refer to equation (3)), the 
exploitation of these half-feasible crossovers must be considered when dealing with 
these kinds of problems. This undertaking requires the modification of the search 
procedure so far discussed, to recognise those crossovers for which only one of the 
resultant vectors is feasible. 

This might be actualised by resorting to a half-feasible crossover only in the absence of 
any that are entirely feasible. Such a scheme respects a preference for the entirely 
feasible variety, but this may be at the expense of an early termination; when no 
feasible crossover is found, the search must complete before allowing the acceptance of 
a half-feasible crossover. This method appears most useful, then, when the number of 
entirely feasible crossovers in G>(P) is reasonably high, yet the chance of discovering 
none cannot be ignored either, say when I P \~ -|(l A I +1). However, this region may be 
quite narrow [7]. 
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An alternative approach, most appropriate for the most highly-constrained examples 
(when I P k< f (I A i +1)), is to regard feasible and half-feasible crossovers uniformly, 
accepting either when unearthed with probability 

_ 1 _ 

?F+H “ («-!)• MS a* (x, y) c P)+ Pr(l E hk (x, y)f]P\= l)]' 

This expression is derived using Theorems 2.2 and 2.3, similarly to that for p F above, 
and can be likewise approximated by keeping a count t F+H of the total number of 
entirely feasible and half-feasible crossovers discovered to date. 

With the restriction that only feasible points are permitted to enter the new population, 
one of the solutions offered by the half-feasible crossover must be discarded. Replacing 
the infeasible offspring with either of the two parents is unacceptable, as this results in 
a forfeit of genotypic information; the values common to the discarded parent and 
infeasible child are lost to the algorithm. 

This amounts to a pre-emptive form of selection, at the micro-scale, for a judgment is 
passed on some of the solution fragments involved, and without regard to the specifics 
of the overall selection mechanism being applied. By sustaining both the parent 
individuals and the new solution produced by the half-feasible crossover, the 
genotypic diversity of the population is maintained. Solutions are thereby appraised 
solely by the selection operator, in complete assent with its intended purpose. This 
requires a flexible population, permitted to grow evanescently in size up to half as 
much again, until selection reduces this in creating the next mating pool. 


4. Feasible Mutation 

The combined influence of selection and crossover, together furnishing most of the 
genetic algorithm's information processing efficaciousness, surmounts that of mutation 
(in contrast, the pure evolutionary strategies [l-2,4-5] rely solely or preponderately on 
the mutation operator). However, with the utilisation of selection and crossover alone 
comes the possibility of losing the ability to investigate some potentially beneficial 
slices of the search space. 

A detrimental loss of knowledge might occur, for example, when some substructure of 
the optimum occurs in combination with other, overwhelmingly poor, partial 
solutions. Over-exuberance of the genetic algorithm in its quest may immediately 
eliminate this individual, and the desired substructure with it, from future 
consideration. Indeed, it may occur that no individual exists in the initial population 
that has some particular component value; perhaps this is a value that is required to 
construct an optimal (or a satisfactory near-optimal) solution. 


10 
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The mutation operator, like its biological inspiration, implies a random alteration to 
some vector component (allele), occurring with low probability. Through mutation, an 
individual may become not viable; in biological systems this results in the individual's 
demise. In the confines of the artificial population, to permit the generation of 
infeasible solutions and then eliminate them is a waste of effort, particularly in 
problems where the feasible region and therefore the chances of retaining viability in 
an unrestrained and capricious dissimilation are small. 

Mutation, like the crossover, should be restricted in respect for the constraint system. 
Towards this goal, consider that mutation is a function 

4*: Z" xZx% -» Z" 


defined as 


H-h 1 1 v x k+\ 

The mutation operator may also be viewed as a degeneration of the crossover operator; 
effectively, the mutated individual is crossed with a random vector, where the 
crossover segment is simply the single chosen allele. For this reason, the search for a 
feasible mutation takes a similar form to that for crossover. 


Theorem 4.1 [7]: 
The mutation 


x'='¥(x,e + x k ,k), where l k -x k <e <u k -x k and ee Z, 


is feasible if and only if, for all /: 1 < i < m , 


( 


a Lk = 0v 


a ik > 0 a £ < 


r(x) t 


a i k < 0 a e > 




*i,k 


Applying this result is simple; given the solution vector and index of the component to 
be mutated, the interval of values of £ that would preserve feasibility can be 
determined, by considering in turn each element of the k th column of the constraint 
matrix A. A negative entry a ik defines r(x) l Ia ik as a lower bound on £, and a positive 

entry describes this instead as an upper bound. 

The possible mutation values so determined, the mutation itself can now be 
performed, by substituting the new value x k + e for the old x k . The calling function is 
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also free to decide whether to permit a trivial mutation (when e = 0), or to try another 
component, or a different candidate solution. 


5. Conclusion 

The resolution by an evolutionary algorithm of a mathematical optimisation problem 
that is constrained to discrete values, and by a system of linear inequalities, has been 
disencumbered by the proposal of a set of new genetic operators that ensure the 
continued feasibility of candidate solutions. This is the emanation of theoretical 
investigations into the behaviour of the traditional two-point crossover operator, in the 
context of a set of viable solutions enclosed by some bounding interval. 

Possible crossover operations applied to points of this feasible set are found to fall into 
two relevant categories: those that exclusively manufacture feasible solutions, and 
those that produce one point that is feasible, and one that is not. The proportions of 
these are understood to vary with respect to the character of the feasible region itself. 

Upon this basis, a crossover operator is propounded that restricts itself to producing 
nothing but feasible candidates. However, this may undesirably limit the information¬ 
processing efficacy of the genetic algorithm, particularly when dealing with problems 
that are highly constrained, for which the number of crossovers solely yielding feasible 
solutions is greatly diminished. Techniques that utilise the half-feasible variety have 
been developed to overcome this, and their effective application alludes furthermore to 
a ephemeral increase in population size. Finally, a mutation operator is offered that 
likewise guarantees the feasibility of its result. 
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